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About 90% of BsBs pairs produced at the T(5s) resonance are initially B^B^ pairs which decay 
radiatively to BsBs- This implies that the BsBs pair will then be in an eigenstate of charge 
conjugation (i.e. C — —1) and therefore in an entangled state. This allows for a determination of 
^ AFs/Fs and the CP phase using a number of possible correlations between the decays of the two 

Bs mesons. In particular, we consider the time integrated correlation; the time ordering asymmetry 
and the time ordering-charge asymmetry, which in addition to time ordering distinguishes Bs from 
Bs, for various combinations of final states. With the statistics of about 0(10^ — 10**) T(5s) events 
available at B factories, we find that the time ordering asymmetry between suitably defined hadronic 
and flavor specific (tagging) decays offers a promising method for determining the width difference. 
The corresponding time ordering-charge asymmetry can also bound the mixing phase. Similar 
observables involving exclusive decays are also considered. At the super B factories with 0(50) 
f~i ■ times greater luminosity time ordering and time ordering-charge asymmetries between inclusive and 

CIh' exclusive modes may also provide additional bounds on the phases in those decays. 
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I. INTRODUCTION 

> 

' Resonance production of B- mesons at electron-positron colliders have proven to be an extraordinarily effective tool 
' to study flavor oscillations. The B-factories to date have largely only used the T(4s) resonance where the only neutral 
CSj B-meson produced is the Bd- To produce the more massive Bg meson it is necessary to operate at the T(5s) peak. 
[ Indeed this mode of operation has already been carried out at CLEO[l| and more recently at BELLE 0, Q where 
' lOOfb^^ of data, amounting to about 10'' BgBs pairs has been collected. 
T— I ' In the last few years, much progress has been made in the study of mixing in the Bs system at DO [3| and CDF 
0^ . the mass difference and width differences have been measured and some bounds have been placed on the oscillation 
■ phase [1, 01 ■ Further progress is expected in the future and also at the LHCb experiment when it takes data. 
^ While these experiments offer the advantage of high event rates, they are limited in the number of final states which 
r '"j ■ they can observe. If new physics contributes to the mixing or decay in the Bs system, it could result in different CP 
rS phases in different final states so doing mixing studies in a larger sample of final states is desirable. 

At a B factory the detector environment is generally cleaner than at a hadronic machine so that a larger set of final 
■ " " ' states can be observed, particularly those involving neutrals in their final state. On the other hand the B factory 
has the disadvantage that the Bs mesons are less boosted so that oscillations with frequency Arris might not be 
observable. Finally, the B-factory does have the feature that the mesons are produced in correlated pairs. The focus 
of this paper will be how to exploit this unique feature to obtain information about Bs mixing. 

In this paper we will consider how the quantum correlations between the Bs mesons allow the determination of 
the width difference of the Bg system through the correlation of inclusive final states between the two meson decays. 
Further studies of the correlations with exclusive or semi-exclusive {i.e. an exclusive state with several quantum 
amplitudes or polarizations) final states could further give information about the width difference and the CP phases. 

We will show that luminosity sufficient to produce 10^ — 10* will generate sufficient statistics to carry out such 
studies however systematic errors originating from the single Bs branching ratios would need to be addressed in order 
to obtain precision results. 

Another way to address this problem is to take advantage of the correlated time evolution of the entangled meson 
pair. We find that time ordering asymmetries between different inclusive states are particularly sensitive to ATg. A 
related time ordering-charge asymmetry can also give the tangent of the mixing phase, in spite of the fact that we 
assume the Am^ oscillations are too rapid to be observed in detail at B factories. This capability will help resolve the 
sign ambiguities in the time ordering asymmetry and correlation measurements which only give the cosine of phases. 
This methodology can be extended to correlations between inclusive and semi-exclusive states to obtain separately 
the CP phase for those decays. The key advantage of using these asymmetries is that they are null experiments which 
vanish in the limit of no mixing and therefore are not subject to large systematic errors from inp ut branching ratios. 

Future B factories are being considered with luminosities 50 times the current machines |9l4l l| . If a fraction of this 
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luminosity is on the T(5s) then mixing effects in a larger set of exclusive decays may be probed including final states 
which are more sensitive to new physics. 

This method is complementary to the methods used at hadronic experiments, in particular CDF and DO where the 
oscillations proportional to Arn^ can be resolved. Note in particular the existing results from DO and CDF[4.. .5., .12\. 

In this paper we will generally be considering observables which should be within the capabilities of modern B- 
factories. For instance a B-factory with a luminosities of O(10^'*cm~^s~^) operating at the T(5s) peak produces 
bb final states with a cross section[2] of (0.302 ± 0.014)nb. Thus, running 1 year there will be O(IO^) such events. 
Upgrade plans for KEK have as the design luminosity 50 ab~^. If about 10% of this luminosity is delivered at the 
T(5s) then this will produce 1.5 x 10^ bb events. 

There are several compelling motivations to preform these measurements. First of all, recent theoretical progress 
suggests that a fairly robust Standard Model prediction of the width difference may be possible and so 
measurement of this quantity could become a good test of the Standard Model. The Standard Model also predicts 
that the mixing and decay phases in the Bg system are small, and the same for all final states with quark content cess 
since the overall phase of this combination is arg{—{VtbV*g)/{VcbV*g)) « 0.02. Thus the observation of a significant 
phase and/or variation in the phase between different cess final states would indicate new physics. 

In Section [n] we discuss the The BgBs correlated state, oscillation and CP violation in decays are dealt with in 
Section IIIII Inclusive and exclusive final states are discussed in Section IIVI time independent and time dependent 
effects follow in Sections |V] and |Vl] respectively. A brief conclusion is given in Section IVIII 

II. THE BjBs CORRELATED STATE 

At the T(5s) peak there is sufficient energy to produce the final states BgBg, B*Bs, BgB^ and B*Bg. The vector 
state decays radiatively, B* — >• jB^, so in all of these cases the final state consists of a BgBg pair and n = 0, 1 or 
2 photons. This means that the charge conjugation of the final BgBg system is directly related to the number of 
photons so radiated[17]. In particular, the initial T(5s) state is in a C = —1 state so if it goes directly to BgBg then 
the meson pair must also be in a C = —1 state. If the transition is through B*Bs or B^Bg then the final BgBg 
pair must be in a C = +1 state because of the associated photon. Likewise if the transition is through B*B* then 
the final state consists of two photons and a BgBg pair so the BgBs pair is in a C = —1 state. Note that this same 
argument applies if the meson pair is produced through a virtual photon (indeed one cannot a priori separate the 
two channels). Current results'2| indicate that about 90% of BgBg pair production is through the B*B* (while the 
relative contributions of the other two modes are less well measured) so the final meson pair is in the C= —1 state at 
least 90% of the time. 

If we let k be the 3-momentum of the Bs in the center of mass of the B^Bg pair, then the wave function is constrained 
by the fact that for a scalar anti-scalar pair P=C=(— 1)^ so that for the odd L, C=— 1 case: 

vI/^=-i (X i= {\Bg{k))\Bg{-k)) - \Bg{k))\Bg{-k))) (1) 
while for the even L, C=+l case the wave function is: 

vI/^=+i (X -L (\Bg{k))\Bg{-k)) + \Bg{k))\Bg{-k))) (2) 

Let us denote a = a{e'^e~ — >■ B^^^ B^^ '')/(T(e+e~ — >■ bb) at the T(5s) peak and let r, r* and r** refer to the fraction 
of this branching ratio which contains 0, 1 or 2 vectors respectively. Current experimental results (T], [3, [3 
a = 19.7 ± 2.9% and r** = 90.114'q%. r and r* are as yet unmeasured but are constrained by the r** result since 
r + r* + r** = 1. Clearly then the BgBg pair is dominantly in the 'Jf^=~^ state, at least 90% of the time. 

The BgBs pairs thus produced at the T(5s) in the C=-l state are therefore similar to the B'^B pairs produced 
at the T(4s). However, because of the different regime of mixing parameters, the quantities which can be measured 
using this effect are somewhat different. 

III. OSCILLATION AND CP VIOLATION IN B, DECAY 

Let us now turn our attention to the time evolution of the Bg mesons. Unlike in the Bd case, direct measurement 
of the oscillations driven by the mass difference are probably not practical at a B-factory since Anis = 17.77 ± 



3 



0.1 ± 0.07ps ^ is so large, it gives rapid oscillations which are hard to resolve. In this paper we will instead focus on 
results related to the component of mixing driven by the width difference. Currently the width difference is not well 
measured, results from CDFi] and D0[1 give \l§\ \ATs\ = 0.154;° °^^ps-^ This is about 20% of the decay rate 
where = l/Fg = 1.472^q q26P^- shall show, this range of mixing parameters allows time integrated and time 

dependent studies to provide information concerning the width difference and CP violation in the Bg system. 

Let us first review the standard formalism for time evolution in a single neutral mesonfio'] and then generalize this 
to an entangled pair of mesons. We denote a state \1/ — a^jS^) +a^i,\Bs) which is a mixture of Bg and Bg by: 



1- 



flip 



(3) 



In this basis we can write the the general mass matrix subject to CPT constraints: 



M 



A 



A 



(4) 



Here A, q and p are general complex numbers where we take the convention that Re{pq) > 0. This matrix then has 
complex eigenvalues fii = A—pq and /i2 = A+pq. We will write these eigenvalues as fii — mi — ^Ti and fi2 = ^2 — -^^2 
where m — ^(mi + and Tg = ^(ri + r2)- The eigenfunctions corresponding to these two eigenvalues are: 



\^i)^\Bi) 



1 



\B. 



1 



(5) 



If ^'(i) is the state at time i, this is related to the state at time=0 by: 



(6) 



where the time evolution matrix U satisfies 



i-U{t) = MU{t) 



(7) 



The solution to this equation is: 



U 



9+ 



9a 



9- 



(8) 



where e5''+*"^ = q/p and g± = | (e"'''!* ± 6"'^^*). 

In the Bg system \q/p\ ~ 1, the current experimental value is[18] \q/p\ — 1.0019 ±0.0047. This deviation from j] — 
is undetectably small for the methods we will discuss in this paper so we will proceed taking the approximation that 
rj ~ 0. This experimental value of \q/p\ is obtained from measurement of the semi-leptonic asymmetry. 



AsLiBg^) 



N{Bg^d ^ t+J^iX) + N{Bg^d ^ i-viX) 



(9) 



using the relation 



62' 



\<ilp\ 



\~ A 

TVa 



1/4 



(10) 



The average value of Ag^iBg) used in obtaining this value of \qlp\ is AshiBg^ = — 0.0037 ± 0.0094 where the authors 
of [13 have combined T(4s) data giving Ash{Bd) with Tevatron data which gives a linear combination of Agj^iB^ 
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with Asl{Bs)- a more resent DO result t20j] which is not included in this average gives the combined asymmetry as 
= -0.00957 ± 0.00251(siai) ± 0.00146(s?/si) where 



A 



SL 



N{b ^ e+i^iX) ~ N{b -> e-veX) 
N{b £+viX) + N{b ^ l-VtX) ' 



(11) 



This is a linear combination of AsL{Bd) and Asl{Bs) which has a 3.2 sigma discrepancy from the Standard Model 
prediction for this quantity: — — 2.3lg'g x 10^^ (predicted). This deviation highlights another utility of T(5s) 
B-factories which should be able to directly measure A\,^ as well as the Bs and Bd components separately and thus 
clarify the comparison of theory to experiment in Asl- 

Let /,; be a single quantum state with Ai and Ai being the decay amplitudes of Bs and Bg to fi respectively. 
Denoting A; — [A^ , Ai\ we will normalize the units of amplitude so that the decay rate to fi for a given initial state 
*(0) is; 



We can rewrite this as 



T,{t)^Tr [uHt) R, U{t) po] 



(12) 



(13) 



where po = ^(0) ^'(0)^ and Ri — A^^A^ is the decay density matrix for a single quantum state. 

Consider now, more generally a state F consisting of several individual quantum states: F = {fi}- The decay rate 
as a function of time thus becomes: 



TF{t)^Tr [uHt) Rf U{t) po] 
where Rp = Ri. In general Rp is Hermitian and thus can be written in the form: 



(14) 



Rf 



Up + vp wpe 
wpe^^^'' Up — Vp 



(15) 



where up, vp and wp are positive real numbers and uj^ > Vp + Wp. If F consists of only a single quantum state then 



ul' — Vp + w"^ 



p — Up -r Wp. 

Let us now expand Eqn. (|14|) for the initial states i?s, B^ and {Bg/Bg} which is an incoherent mixture of an equal 
number of B^ and Bs- The corresponding density matrices for these initial states are: 



PiBs 



1 





PiBs 




1 



pi{Bs,Bs}) 



1 
1 



The time dependent decay rates of these initial states to F in the limit 7/^0 (i.e. \q/p\ — 1) are: 



where 



and 



T{Bs ^ F) 
T(Bs ^ F) 
T{{Bs/Bs} F) 



e ' [upChy 



e [upChy 
[upCh 



V 



- vpC'x - Wp {Shy cos{(f) + 9p) - Sx sin(0 + 6p))] 
vfCx — Wp {Shy cos(0 + 6p) + Sx sin(0 + 9p))] 
wpShy cos{4> + Op)] 



Xs = ATOs/Fs 



Ar,/(2r, 



(16) 



(17) 



(18) 



Cx = cos{xsVst) Sx = s\i\{xsTst) Chy = cos\v{ysVst) Shy = smh{ysTst) 



(19) 
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The branching ratio to a particular final state is the time integral of the above. In particular we denote the branching 
ratio from the initial state {Bs/Bg} by 

~ UF- WFVs cos(0 + 9) 
= O^^s 

which is the average of the branching ratios of Bi and B2 to this final state. 

Let us now extend the above formalism to the system of a correlated BgBs (i.e. produced at the T(5s)) state where 
the meson 1 with momentum k decays to state Fi and meson 2 with momentum —k decays to final state F2. We can 
write time dependent decay rate as: 

r^,F,{tl,t2) = SF.F.Tr [{U\h)RF,U{h))Z±{U\t2)RF,U{t2)fZl\ (21) 

where the superscript T stands for transpose, SfiF2 is a combinatorial factor; SfiF2 = 1 li Fi ^ F2 and SfiF2 = 5 if 
Fi = F2 and 



1 
±1 



(22) 



is the matrix representation of the initial wave function, ^''^=='=1. Here ti is the time of decay for meson #1 and t2 is 
the time of decay for meson #2. 

Expanding the above for a C = — 1 initial state, the result is: 

r^;^p.,(ii,t2) = 25i.,F,e-r-(*^+*^) [{uiU2 - tmCtct)C^y + {uiW2Ct - U2WiCt)S^y 

+ {v2WiSf - ViW2S2)S~ - {viV2 + Wl Ui25'f )C~] (23) 

while for an initial C = +1 state, the decay rate is: 

r+ p,^(tl,i2) = 2SF,F2e-^'^*'+*'^ [{utU2 + W^W2CtCt)C+y - {U2WiCt + UiW2Ct)S+y 

+ {V2WiSt + ViW2St)S+ - {VIV2 - WiW2StSt)C+] (24) 

where 

Cf = cos(</) + 9,) Sf = sin(<;A + 9,) 
Ct = cos((ti ± t2)x,T,) = sin((ti ± t2)x,T,) 

C± = cosh((ti ± i2)2/«r,) 5± = sinh((ti ± t2)2/«r,) (25) 

Integrating these results over ti and t2 we obtain the correlated branching ratios 



B-{F,F2) = ^ 



(UlU2 - WiW2CfC^) ^ - {viV2 + WiW2SfS^) } 

B+{F^F2) = [{u^U2 + w,W2Cfct)j^^^ - {u2W,Ct + u,W2Ct)j^^,^ 



(26) 
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+ {v2W,St + V,W2St)^^f^, + {V,V2 - W,W2StSt)^:^] (27) 

where these are the ratios with respect to the total number of B^Bg pairs produced in each of the two CP states. The 
fraction of bb events at the T(5s) peak will therefore be: 

BUF1F2) = a ((1 - r*)B-{F^F2) + r*B+{F,F2)) (28) 

where we assume that no attempt is made to distinguish between the C = +1 and C = — 1 BgBg pairs. If these 

cases can be distinguished, for instance by counting the number of photons associated with the system, then an 
improvement in the statistics may be obtained although such a potential improvement is somewhat limited by the 
fact that at least 90% of the meson pairs are in the C = — 1 state. 
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IV. INCLUSIVE AND EXCLUSIVE FINAL STATES 

The key to obtaining basic physics parameters from the correlations and asymmetries we will discuss below is to 
choose final states where there is some a priori knowledge of the mixing strength wp defined in Eqn. (jlSp . To this 
end we will consider two opposite limits in which that is the case. On the one hand we will consider inclusive states 
which means a large fraction of the decay modes. Depending on how you select such modes, wp can be either 
or related directly to y. On the other extreme CP eigenstates or related exclusive states which provide a case where 
\wf/uf\ = 1- 

Here we will consider two categories of inclusive final states, first of all flavor specific "taggable" final states which 
exclude quark content cess. Second of all "hadronic" final states which include cess. 

The set of taggable decays includes all decays where the flavor of the meson can be determined from the decay 
products. For example in the decay Bg — ?> fi~^i'D~ it is known that at the instant of decay that there was a Bg (i.e. 
specifically a bs state) and not a Bg (i.e. bs). We will denote taggable decays which indicate an initial state of Bg 
by t+ and taggable decays which indicate an initial state of Bg by t—. If we are not concerned with the flavor of the 
initial state we will denote the state as t = {t+) U (i— ). For such decays wt = since regardless of whether they tag 
B = +1 or _B = — 1 they cannot mix between Bg and Bg. 

These modes are a signiflcant fraction of Bg decays. They consist of all semileptonic decays as well as most hadronic 
decays which do not have quark content cess. For instance we can include many hadronic decays containing only one 
charmed meson. It is advantageous to be able to include as many decays in the taggable sample as possible; we will 
assume that the tagging rate is 30% of all Bg decays. 

If there is a partial rate asymmetry in taggable decays then vt ^ 0. Initially we will not generally be considering 
observables which are particularly sensitive to this kind of effect. 

The category of "hadronic" decays, which we denote "/i" may include all hadronic decays or, more generally a 
subset of hadronic decays that includes decays with quark content ecs's. To optimize the utility of this sample, it is 
best to include all cess final states and as few other hadronic states as possible into the "hadronic" sample. As we will 
show below, cuts which are tight enough to reduce the hadronic sample by about 20% while passing all cess states 
greatly improve statistics in some cases. 

For hadronic states Vh would indicate that there is a partial rate asymmetry. As with the taggable decays, the 
observables we discuss in this paper will generally not be sensitive to so we will take Vh = 0. Wh however will be 
non-zero and is, in fact, tied to the Bg lifetime difference since the lifetime difference arises from a difference in the 
decay rate of the eigenstates to ecs's states. Note that this also leads to a phase 6h. 

Exclusive flavor neutral states which consist of a single quantum amplitude, such as -F = DfDj, allow the direct 
measurement of the quantity ycos{4> + Oi). For such states wp — y^uj^ ~ Vp so assuming that the partial rate 
asymmetry is not large, wp ~ up is a good approximation. Thus, the matrix Rp depends only on up and the CP 
phase 9p. In the limit of CP conservation this is further constrained. If F is a CP eigenstate and CP = +1, up = wp 
and 9p — vp = while if _F is a CP = —1 state, up = wp, 6p = tt and vp — 0. 

Other exclusive states such as F = D'^*Dj* consist of multiple quantum states (in this case due to polarization). 
The parameter Wi is therefore not constrained. We may however be able to obtain information about the relative 
contribution of different amplitudes that make up F through the study of {Bg/Bg^ F. For instance in the case 
of F = Df*Dj* we can learn the contributions of the different polarization states through studies at a hadronic Bg 
experiment. 

More generally, a semi-exclusive state consists of a small set of inclusive states such as F — tp + X. Most likely 
there is no simple way to determine wp or 9p. However, if you combine information from time ordering asymmetries 
and time ordering-charge asymmetries discussed below, you can obtain the phase of such a decay. The Standard 
Model implies that all such phases will be small so if a large phase is discovered in any semi-inclusive set, this could 
be evidence for new physics. 

V. TIME INDEPENDENT CORRELATIONS 

Let us first consider the effect that mixing has on the time independent correlations between final states. If no 
mixing were present, the null hypothesis, the overall branching ratio to the state FiFj would be given in a simple way 
by the product of the branching ratios for each of the Bg meson giving 

B'^fiFF,) = 2aSp^p^B{F,)B{F,) (29) 

With mixing present, we will have a deviation of the measured value of B^g from this expectation. The magnitude 
of this deviation will thus tell us about the mixing process. 
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In order to carry out this program, however you need to have an accurate value for the basic nuU hypothesis so 
B{Fi) and B{Fj) must be well determined. The systematic error in B2f\F,Fj) is therefore likely to be the mam 
limitation in using this technique to probe mixing. 

Let us first consider the application to inclusive states and so apply Eqns. (^01 [?fl and 1^51) to t and h states. 

In the case of taggable states, since Wt = 0, 



The origin of AFg is the rate difference within the hadronic decays so that y is related to the hadronic decay mixing 
parameters by: 

ysTs = WhC0fi{(j) + 6h) (31) 

it follows then that 

~ Uh~whyscos{(l) + 0h) 

= ^li^) ^''^ 

It is convenient to write the correlated branching ratios strictly in terms of B since these are separately determined 
experimental quantities. So turning now to the correlated branching ratios for tt, th and hh final states we obtain. 

B5.(tt) = aBUl - (1 ~ 2r*)yj) + 0{x-') 
B,,{hh) = aBl{l-{l-2r*)yl{B-'~lf)) + 0{x;') 

BUth) = 2aAs,,(l + (l-2r*)y2(B-i-l)) + 0(a:;2) (33) 

Here we drop the 0{x~^) terms since Xs is large for the Bs system. 

In each of these cases, the correlated branching ratio depends only on y1 and other measurable branching ratios 
and is independent of the CP violating phases. In TableUwe show the number of bb events required at the T(5s) peak 
to give a 5 — cr statistical deviation from y = both in the case of ys = 0.1 and ys = 0.05. For the hadronic states, we 
also consider the scenarios where B/j = 0.7 and B^ = 0.5. From these results we see that applying cuts which reduce 
Bfi from 0.7 to 0.5 will be very helpful in determining y^ with this strategy. For instance in the hh case, A^5s(5(t) is 
lowered from 120 x 10^ to 8 x 10^. Note also that since hh, ht and tt correlations all measure y^ it makes sense to 
combine the results from each of these combinations of final states. Combining data in this way will also lead to a 
reduction in A^5s(5(t). 

Thus there may well be adequate statistics to carry out this program, either presently or in the near future especially 
since BELLE as mentioned previously, already has accumulated appreciable amount of data at the T(5s). However 
the fractional deviation of the B^s from -B™" is 0{y^) ~ 1%. This implies that the accuracy of the input values for 
Bt and Bh needs to be less than 1% for the signal to be observable. 

In the limit of CP conservation where the states are CP eigenstates, Eqns. l\'6'S\i can be understood in terms of a 
simple argument. Consider the case of B^sitt). The decay rate of each of the eigenstates to taggable final states is 
the same so the branching ratio will be inversely proportional to the total decay rate. Thus 

B{B,^t) = —^^ B{B2^t)= (34) 

rs(l-?/s) Tsil + ys) 



so taking the average, we obtain Eqn. (|30l) 

If the initial state is C = — 1, then the two meson state is (|_Bi)|i?2) — \Bi)\B2))/V2 so that B^ {tt) — B{Bi — > 
t)B{B2 t). Likewise, the C = +1 state is {\Bi)\Bi) + \B2}\B2})/V^ so B+{tt) = \{B{Bi t)B{Bi t) + B{B2 
t)B{B2 t)). Putting in the eigenstate branching ratios Eqn. (p4|) into these expressions we can thus derive the tt 
correlation in Eqn. (1551) . 

To obtain the other expressions in this limit, the hadronic branching ratio of the two CP eigenstates is 



(35) 
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which leads to the other two correlations in Eqn. (155)) . 

Let us now consider the correlations involving an exclusive or semi-exclusive final state Y. For such a final state, 
let us define the quantity 

WY cos((/» + 6ly) 

Py = (36) 

Uy 

which is the quantity that we wish to measure. It greatly simplifies the expressions below to define the related 
quantity: 

" - ^^-^^ (37) 



1 - PyVs 

In the limit that CP is conserved, Py is 

Py = (38) 

where Y'^ is the subset of Y which is CP = +1 and Y~ is the subset of Y which is CP = — 1. More generally, for 
taggable states, Pt = while for hadronic states Ph — ysB^^ /{I ~ y^); thus: 

Pt = -y. Ph = {Bi^^ - (39) 
Using this notation for the branching ratio to the state Y: 



By = V^f" (40) 



The correlation between Y and the h and t inclusive states are: 

B^ty) - ^aBtBY (l + (1 - 2r*)y,Py) + 0{x-^) (41) 
B^sihY) - 2aBHBY (l - (1 - 2r*)y,(B,;i - l)Py) + 0{x-^) (42) 
The correlation between two different states in general is: 



BUy^yJ) = 2aSYY.,BYBY, V ' ' ^^IPy^Py,) + Oixj') (43) 

As discussed injl^, [2^, it is easy to understand these correlations in the limit of CP conservation. In particular, 
suppose that F is a CP=-|-1 eigenstate so Py — +1- If we start with an initial charge conjugation —1 BgBg state, 
then if one of the mesons decays to Y, the other must therefore be in the B2 state. The probability of it decaying to 
a taggable state is therefore B{B2 t) given in Eqn. Conversely if we start with an initial charge conjugation 

+1 BsBg state, then if one of the mesons decays to Y, the other must therefore be in the Bi state. In this case 
the probability of decaying to a taggable decay is B{Bi — )■ t). From this we can derive the tY correlation above in 
the case Py = 1- We can generalize this argument to a case where \Py\ < 1 considering separately the Y~^ and Y^ 
components. 

In summary then the correlations between inclusive states in Eqn. (I33p can determine \ys\ while the correlations 
between Y and h or t give Py which in turn gives ywY cos{(f) + Oy)- If wy can be determined in some way, then the 
cosine of the phase is determined. 

If Y is flavor neutral exclusive state (i.e. a CP eigenstate), then wy — uy- For states with multiple amplitudes 
such as D'*D*, -00, DsDsT]', w must be determined from detailed analysis of the final state or the state needs to be 
separated into its constituent quantum states, for instance by polarization or Dalitz plot analysis. For semi-exclusive 
states such as (/) -f X there is no way with this kind of data to factor Py into wy and cos(0y + </<), so more information 
is required to do this. 



9 



Final State 


Inputs 




JV5s i^ocr j/iu \{j — v.yjo) 


T(5s 


— > 




Bt = 0.3, r* = 


0.1 






91 


OOU 


T{5s) 


ht 




Bh = 0.7, Bt = 0.3, 


r* = 


0.1 




25 


410 


T(5s) 


ht 




Bh = 0.5, Bt = 0.3, 


r* = 


0.1 




6.5 


100 


T(5s) 


— ^ /i/i 




Bh = 0.7, r* = 


0.1 






120 


1900 


T(5s) 


— hh 




Bh = 0.5, r* = 


0.1 






7.8 


130 


T(5s) ^ 


hh; ht; tt 




Bh = 0.7, r* = 


0.1 






16 


260 


T(5s) ^ 


hh; ht; tt 




Bh = 0.5, r* = 


0.1 






4.9 


78 


T(5s) 


Yt 


Bt 


= 0.3, By = 10-^ Py 


= 1 


r* 


= 0.1 


33 


130 


T(5s) 


^ Yh 


Bh 


= 0.7, By = 10"^ Py 


= 1 




= 0.1 


76 


300 


T(5s) 


Yh 


Bh 


= 0.5, By = 10~^ Py 


= 1 


* 


= 0.1 


20 


78 



TABLE I: The value of N5s{5<j), the number of bb events at the T(5s) required to observe a 5cr deviation from the null 
hypothesis Eon. 1291 assuming perfect knowledge of the input branching ratios, for the pairs of final states indicated. The results 
are shown for j/s = 0.1 and ys = 0.05. 



Final State 


Bi/Bj 


a,/a,(%) 


iV5.(5a) (10») 


DtD^/DtDj 
K+K-/K+K- 
DtD-/K+K- 

7V°KsKs/7V°KsKs 


1.1 X lO^Vl.l X 10"^ 
3.3 X 10-V3.3 X 10-"' 
1.1 X 10"V3.3 X 10"^ 

9 X lO-'/g X 10"'' 
3.7 X 10"V3.7 X 10"^ 
9 X 10"''/3.7 X lO"'^ 


1/1 
100/100 
1/100 

4/4 
25/25 

4/25 


12 
130 
20 
110 
1600 
210 



TABLE IL The number of 66 events (unlike the previous Table, here in units of 10^) required to distinguish between the 
maximum and minimum possible correlation between two CP eigenstates at the 5(j level including the acceptances for the final 
states shown. 



Correlations between two semi-exclusive states may also be used to determine Py by using the correlation in 
Eqn. (USD. 

In the case where the two states are the same or have the same CP eigenvalue this method has the advantage that 
it is almost a null experiment and it is therefore not subject to the contamination due to systematic errors in the 
input branching ratios. Looking at Eqn. (^5)) we see that in the limit of r* — )■ and Py^ — Py^ = ±1 which would be 
the case if CP were conserved then would be 0. This limit can be understood in terms of Bose statistics. In the 
C — —1 state, a Bg pair consists of one CP=-|-1 (i.e. Bi) and one CP=-1 (i.e. B2) meson so you would never see two 
decays to the same CP eigenstate. 

This correlation is not directly sensitive to y but it is sensitive to the mixing angle of the two states. As an 
illustration, let us consider the case where we have two states with B = 10^^ where r* = 0.1 that are CP=+1 
eigenstates. If there is no CP violation so Py^ — Py. = +1 then i35s(no mixing) = 8 x 10"*; where the fact that it is 
non-zero is due to the term proportional to r* . On the other hand, suppose that there were large mixings in one or 
both of the channels so that \Py.Pyj \ << 1 and so i35s(large mixing) = 4 x 10"^ about 5 times larger. Thus if you 
have TV = 7.5 x lO*" events you can rule out the large mixing scenario at 5a. 

Most likely, however, one can probably not find a CP eigenstate decay mode with a branching ratio this high if the 
acceptance is factored in. In Table |lT] we consider the number of bb required to distinguish between the minimum and 
maximum possible correlations at the 5a level. The CP=-|-1 states we consider in particular are DfDj where we 
assume that the acceptance of this final state is a = 1%, and K~^K~ which should have acceptance nearly a = 100%. 
The latter state has the advantage that it has a significant penguin contribution and so is more likely to be influenced 
by new physics. The CP=-1 states we consider in particular are iprj where we assume that the acceptance of this final 
state is a = 4%, and Tr^KgKg where we assume the acceptance is a = 25%. 

We can see that the numbers are large even for super B-factories. Perhaps the most promising cases are 
DfD^/DfD~ and DfD~ /K~^K~ where the limiting factor is the acceptance of the DfD~. 

VI. TIME DEPENDENT EFFECTS 

Looking at Eqns. (^51 [M)) a couple of features of the time dependence are apparent. Let us denote t± = ii ± ^2 so 
that the expression for r^{ti,t2) is of the form 6"^'"''*+ K we integrate this over t+ we obtain 
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where t_ ranges from ~oo to +cx) and the superscript (— ) on B indicates a C-odd initial state. 

The fact that t+ integrates out triviahy is, of course, exploited in the design of the _B-factory. We can determine 
the function / by binning events according to t- and taking into account the exponential prefactor. Because the B 
mesons are created with proper motion in the lab frame, i_ can be inferred by the physical separation between the two 
decay vertices. In particular, the original interaction vertex where the mesons are created need not be determined. 
This is a useful feature of F" since the e+e" interaction vertex cannot be directly observed. 

In contrast, the expression for r+(ti, ^2) is of the form 6'^"*+ g{t+). If we integrate out the variable t_ we obtain 



= t+e-^^^'^git+). (45) 

where t+ ranges from to +00. To determine t+ we do indeed need to know the location of the interaction vertex. It 
is thus more difficult to study the time dependence of . In fact experimental studies of the time dependence will 
be greatly helped by the feature that the meson pair is in a C = —1 more than 90% of the time because the time 
difference is easier to measure than the time sum (t^) at an asymmetric B factory. 

The terms proportional to 5*,^^^ and S~ in Eqn. (j23p have the property that they are antisymmetric under ti 0^2- 
An observable which has the same symmetry will therefore be sensitive to these terms. The — sm{Ams{ti — ^2)) 
term oscillates at the rate Am^ so it is not readily observable at B factories by directly resolving the oscillations. 
However, the asymmetry A[j discussed below does offer the prospect of sensitivity to this kind of time dependence. 

The simplest such observable is the time ordering asymmetry: 



" \ -1 if < h ^ ^ 

the expectation value of A^j will receive contributions from the coefficients of S^y and which we denote: 



Ql = 2iu,w,Cj - u.w.Cf) = 2TiB,B,{l - y'){P, - P.,) Q^^ = 2{v,w,Ct - v,w,Cf) (47) 

where for flavor neutral final states with Vi = Vj = 0, the term proportional to Qf^ vanishes. 

The expectation value of Aij is, of course, just the asymmetry between the two decay modes according to which 
decays first. Since the two meson final state at B factories has a proper motion in the lab, this asymmetry can be 
calculated by: 

(cases where F; happens later) — (cases where Fj happens later) , , 

< > - all F.F, events ' ^''^ 

where "happens later" translates into "decays further downstream" . 
Evaluating this expectation value we obtain in the fiavor neutral case: 

< A, >= (1 - r*) (qIj^) = (1 - mP, P)ys (49) 
The number of bb events at the T(55) required to observe this asymmetry with a n — a significance is 

Tl^ 1 Tl^ 

A^ 2aBiBj 2a(l - r*)^B^Bj{Pj - Pify^ 

Let us now specialize this formula to the inclusive states. Since the two final states must be different to form this 
asymmetry, for a th final state this becomes: 
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NUna;A,,{th)) = — -^-^ (51) 

2a(l - r*YBtyi 



where the dependence on Bh in the numerator comes from the expression for Ph in Eqn. [ 
Likewise, if we correlate semi-inclusive state Y with i or ft, states: 



N5s{na]tY) ^ ^ . (52) 

2a{l~r*yBtBY{PY+ysYvl 

9 

n 

Nzs(n<j;hY) = . ^- (53) 

2a{l~r*fB^By{PY+Vs{Bi,'-l)?yl 

In Table Hill we show the results for A^5s(5ct). Again we consider a generic state Y with branching ratio 10~^ and 
Py near 1. 

Looking at the value of A^5s(5cr) for the asymmetry in the ht final state, we see that the asymmetry Aij requires 
lower statistics to measure y^. This is largely due to the fact that the taggable decay is more likely for the long lived 
B2 state while the hadronic decay is more likely for the short lived Bi state and the two tendencies combine in the 
asymmetry since the dominant C — ~1 state is always B1B2 by Bose statistics. 

Furthermore, this is a null experiment, absent mixing the asymmetry will be and so there is no large systematic 
error brought in due to the uncertainty of the input branching ratios. 

Using the method of [l^l one can devise an observable which measures this term with optimal statistical efficiency. 
If we neglect terms of 0{r*) and 0{xj^) then the optimal observable to measure the 5*^^^ term is 

T- = tanh(ii - t2) (54) 

This time dependence is proportional to the ratio between the time ordering asymmetric term in Eqn. (|23|) (oc Shy) 
and the time ordering symmetric term (cx Chy)- For small ys this observable offers some improvement, about factor 
of 0(2), over the unweighted time order asymmetry. 

Consider now the case where taggable states are correlated with some other flavor neutral state, either h or Y . In 
this case we can take into account the sign of the tag and define the CP odd time ordering-charge asymmetry 

A't, = BAt, (55) 

where B is ±1 for t± taggable states, i.e. events where the tagging decay indicates a Bg are weighted +1 while the 
events where the tagging decay indicates a Bg are weighted —1. 

Since Vt+ = ut+ = Ut- = —Vt-, this asymmetry is sensitive to Q^j. Looking at the definitions of QJj- and Q^j we 
see that 

Qf,±)^. = ±tan(0 + 0,)Q^,^)^. (56) 

therefore 

<At, > = il-r*){P,+ys)yg (57) 
< A;^- > = (1 - r*){P, + y,)il_^ tan(0 + 6,) (58) 

L + Xg 

If both of these asymmetries are measured then tan(0 + 9j) can be determined from the ratio without having to 
separately measure wj. In practice the Standard Model prediction for tan(6' + </))« .02 so this method will generally 
bound (or discover) large phases due to the presence of new physics. 

There are three distinct ways in which this pair of asymmetries {A and A') may be used. 

1. With the inclusive states, Ath can be used to determine \ys\. A[f^ can be used to find tan(0 + 9h). 
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2. With an exclusive state Y, Aty determines cos(0 + 9y) (see Eqns. [551 and [57| and then A'^y separately gives 
tan(0 + 9y). 

3. For a semi-exclusive state such as Y = cf) + X or Y — ijj + X obtain tan((/) + from the ratio of Aty and Afy 

In all cases of phase measurement, a large phase or discrepancies in phase between different modes indicates new 
physics. 

In Table IIIII we show some sample calculations of the asymmetries that might be seen in various decay modes and 
the statistics required to obtain a 5cr signal for the asymmetry. For each combination of modes we use an assumed 
acceptance a and give product branching ratio (factoring in a) 2aaBiBj with respect to the total number of bb events 
at the T(5s) peak. We then give an estimated time ordering asymmetry Aij which allows us to calculate the number 
of events required for a 5<j signal. 

For combinations with taggable decays, we can use Eqn. (1581) and estimate the time ordering-charge asymmetry 
A'^j assuming tan(0 + 6i) = 1 and the corresponding A^5s(5(t). 

We can compare these Nc^s{^<^) values for the two asymmetries with the ^ 1 x lO'* which could be typical of current 
B-factories and 1.5 x 10^ for a 5 ab~^ super B-factory (i.e. assuming a 50 ab^^ luminosity with about 10% of the 
running devoted to the T(5s)). 

For the case of the inclusive combination ht both asymmetries may be within the range of current B-factories. We 
have also included combinations of inclusive modes with exclusive modes which would likely require a super B factory. 

First of all there is the exclusive mode DgDs which, in the standard model, should be sensitive to the same phase as 
the ht combination. We have considered it with a 10% and a 1% acceptance where in the latter case somewhat more 
than 10* bb events are required. The K'^K~ mode considered would have both tree and penguin contributions hence 
new physics in a QCD penguin could contribute there. The case of ilxj) has already been studied through oscillations 
at DO and CDF. Polarization analysis is helpful in separating the CP even from CP odd amplitudes. In the cases 
where we consider inclusive and exclusive modes with {Ks'k^)d: i-e. a Dq which specifically decays to Ks'k^ , the final 

state connects and D so in the Standard Model there is sensitivity to the CKM phase 7. The more inclusive 
states tp + X and (j> + X have larger branching ratios and the phase in those modes should agree with the overall 
mixing within the SM. A discrepancy of such phases would therefore indicate New Physics. 

VII. CONCLUSION 

In conclusion, with a sample of 0(10"^ — 10*) bb events at the T(5s) peak there is the prospect of making a precision 
determination of AFs/Fs through the study of tt, hh and th correlations. This will allow for the testing of the 
Standard Model prediction of the width difference. 

Time independent correlations between various combinations taggable and hadronic decays have the disadvantage 
that there is a large systematic error originating from the input branching ratios. This can be remedied by using 
time dependent observables. One promising observable to use is the time ordering asymmetry between hadronic and 
taggable decays. The time ordering-charge asymmetry, which in addition to time ordering requires distinguishing Bs 
from Bs, also can constrain the mixing phase although much more statistics would be needed to measure the expected 
Standard Model value. 

At super B factories with about fifty times more luminosity, it becomes feasible to consider time ordering and time 
ordering-charge asymmetries with exclusive states and taggable or hadronic decays. Choosing specific exclusive decay 
modes can thus target different physics issues. 

Time independent correlations between two exclusive CP eigenstates are not subject to the large systematic errors. 
For branching ratios about 10^"^ these correlations can be sensitive to large CP-phases with B factory statistics. 
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Final State (Fi/Fj) 


a=Acceptance(%) 


2aaBiBj (10"'') 


A,, (%) 


iV5s(5a;A) (10'^) 


(tan(,^ + e) = l) (%) 


N^si^a-A') (lO'^) 


h/t 
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84000 
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1.8 


0.49 


12.6 


DtD-/t 
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120 (12) 


10 


21 (210) 


3.78 


150 (1500) 


DtD7 Ih 

s s 1 
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280 (28) 


8.6 


12 (120) 






K+K- /t 


100 


4.0 


10 


630 


3.8 


4400 


K+K-/h 


100 


9.2 


8.6 


367 






ipcj)/t 


10 


11 


5.6 


725 


2.1 


5100 




10 


26 


4.8 


417 






{KsTT°)D(t)/t 


10 


0.37 


10 


6800 


3.8 
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{Ksii°)D<j>/h 


10 


0.84 


10 
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^ + x/t 


10 


140 


1 


1800 


0.38 
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ip + X/h 


10 


325 


1 


770 






<t) + X/t 


50 


1700 


1 


150 


0.38 


1000 


<}) + X/h 


50 


3900 


1 


63 






{Ks7v'')D+X/t 


10 


923 


1 


271 


0.38 


1900 


iKsn'>)D + X/h 


10 


2200 


1 


116 








49 


2.7 


10 


37 


3.8 


260 


(j)'y/h 


49 


6.3 


10 


16 







TABLE III: This table shows an estimate of the requirement for observing the time ordering and time ordering-charge 
asymmetry of the given pairs of final states shown in the first column. For each final state pair, we have assumed the 
acceptance a. Note that for the Ds pairs we have taken a 10% and a 1% scenario (shown in brackets). In the third column we 
show the effective branching ratio compared to e^e~ — >■ hh. In the Aij column we show the estimated time ordering asymmetry 
assuming that cos{4> + 0) ~ 1 as expected in the Standard Model. For the CP eigenstate, D^DJ and K^K~ we assume that 
\P\ = 1, i.e. no large direct CP violation in the decay. For the inclusive states + X and (p + X we suppose the value of 
P ~ 10%. The next column shows A'^5s(5ct), the number of bb events needed to observe the given time ordering asymmetry with 
5a statistics. For the cases where a decay is with a tagged decay, we give the time ordering-charge asymmetry Aij assuming 
that the total phase is tan(<^ + 9) = 1 as might be the case in New Physics. Likewise, we also give the number of bb events 
needed to observe this time ordering-charge asymmetry with 5a significance. 
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